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Solitonic solution-generating methods are powerful tools to construct nontrivial black 
hole solutions of the higher-dimensional Einstein equations systematically. In five dimen- 
sions particularly, the solitonic methods can be successfully applied to the construction of 
asymptotically Minkowski spacetimes with multiple horizons. We review the solitonic meth- 
ods applicable to higher-dimensional vacuum spacetimes and present some five-dimensional 
examples derived from the methods. 
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§1. Introduction 

1.1. Outline 

In the legendary era of general relativity from the 1970s to the early 1980s the 
solution-generating techniques have been greatly developed in four dimensions and 
applied to construct new series of axisymmetric stationary solutions extensively: See 
Ref. Hj) for a comprehensive review. Particularly the practical and useful solitonic 
methods based on the integrability of axisymmetric and stationary gravitational field 
equations were derived by several authors. For example, the Backlund transforma- 
tion analogues were invented by Harrisoil^ and Neugebauer,® and also a certain 
kind of inverse scattering method was developed through the successive works by 
Belinsky and Zakharov.SI' 

In short, in both the methods a known solution (i.e., a seed) is transformed 
into a new one by adding 'solitons' which have formally a role analogous to ordi- 
nary solitons appeared in general nonlinear physics. In these methods the solitons 
can be added in principle with arbitrary number so that a sequence of an infinite 
number of solutions is generated systematically. In axisymmetric and stationary 
cases, nontrivial solitons can bring rotations into the gravitational systems. As one 
of remarkable applications of these methods, Kramer and Neugebauer succeeded to 
generate so called double-Kerr solution^ as vacuum solutions of Einstein equations, 
and opened a new way to the study of gravitational systems with multiple horizons 
since the work done by Israel and Khan based on the Weyl solutions.!^ In four 
dimensions the multi-horizon vacuum systems in asymptotically flat spacetimes gen- 
erally have some flaws like conical singularities. So introducing spin-spin repulsion 
effects by adding nontrivial solitons has been expected to remove such flaws. It 
seems, however, that such completely regular multiple-black-hole spacetimes can not 
be realized in four dimensions.''' 

Mainly in the 1980s, several pioneering works of more than four dimensional 
spacetimes have been done. They have also introduced the solitonic methods to 
analyze the highe-dimensional gravity in the context of the standard Kaluza-Klein 
theory with compactified extra dimensions, and further developed studies of higher- 
dimensional gravitational objects like K-K black holes to some extent.^' However, 
the pace of advancing the study was rather slow because of the increasing intricacy 
of the methods in higher dimensions and also the absence of new stimulus after the 
rise of the unified theories in particle physics including the first revolutionary period 
of string theory. 

The recent discovery of the solutions of five-dimensional black ring, especially 
rotating regular black ring^' has injected new life into all the related subjects of 
higher-dimensional gravity. Inspired by this, systematic construction of new solu- 
tions of five-dimensional black holes has also progressed remarkably using solitonic 
methods. 

The recent progress has been done mainly along two ways, though mutual inter- 
actions between them often occurred. First the solitonic method based on the Neuge- 
bauer's Backlund transformation—' was applied to the five-dimensional Mikowski 
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spacetime as an appropriate seed, and the solutions of black ring with rotation in 5^ 
were discovered P^''^^!' This method was successively applied to several seeds, and 
generated several interesting solutions (5^-rotating black ring^^' and multiple black 
hole systems like black di-rin^HI). However, its application is rather restricted due 
to the metric form adopted to make the method available. Actually, as described in 
section 3.1 the method can only deal with the spacetimes whose metrics have only 
one off-diagonal component. That is, the black holes with two independent rotations 
cannot be constructed by this method. Soon after the above attempt the original 
Belinsky and Zakharov's inverse scattering method was also used to generate some 
of the above solutions .'^''^ However, the original inverse scattering method also 
suffered from the same restriction as that of the former method (see §2.1). 

Following the above works an ingenious device was introduced by Pomeranskjff^ 
to make the inverse scattering method fully effective in higher dimensions ( see §2.2 
for details of this idea). In this work five-dimensional Myers and Perry black holes 
with two rotations were reconstructed. Once the powerful method given, several 
interesting solutions have been generated with the inverse scattering method. The 
multiple black hole systems with single rotation, so called black Saturn and black 
di-ring, have been also constructed by this method, '^''^ though these solutions can 
be derived using the former methods because the metric sector corresponding to the 
three Killing vectors has the structure of 2+1 block matrix. Moreover, the solution of 
a double Myers- Perry black hole is constructed by the inverse scattering method. 
The inverse scattering method shows its power when the off-diagonal components 
generally appear in the metric sector corresponding to the Killing vectors. That is, 
when the doubly spinning rings or rotational black lenses must be considered. The 
solutions of doubly spinning rings and rotational black lenses were constructed in 
Ref. [2T]) and Ref. [22]) . respectively. As another nontrivial example, the orthogonal 
black di-rings (bicycling black rings) were also constructed and the regularities and 
other physical properties were investigated.'^'^' 

The solution-generating tasks will continue, because there remain a lot of in- 
teresting solutions to be dealt with by using the solitonic methods. So at present 
it may be useful to present a summary of solitonic methods and exhibit some of 
the generated solutions. In the following, for short of space, we concentrate on 
solitonic solution-generating methods corresponding to higher-dimensional vacuum 
Einstein equations and the axisymmetric and stationary spacetimes that asymptot- 
ically behave like a Minkowski space. We just comment here that there are several 
works where they have tried to generate the solutions of other systems like Einstein- 
Maxwell ,123 'E3 Einstein-Maxwell-Chern-SimoiPJ and Einstein-Maxwell-DilatoiJ23 

systems. After a brief introduction of the rod structure below, §§2 and 3 are devoted 
to describing two different solitonic methods mainly used so far: the inverse scatter- 
ing method and the Backlund transformation analogue, respectively. In §§4 and 5 
the procedures to generate solutions with single horizon and disconnected horizons 
are demonstrated, respectively. Here, only the procedures of the inverse scattering 
method are shown. 

Notations and Conventions 
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Through this chapter we adopt a higher-dimensional canonical metric form to 
describe the L'-dimensional stationary spacetimes with D — 2 commutable Killing 
vectors (see the metric form (j2-l|) in §2.1). The coordinates associated with the 
D — 2 Killing vectors are denoted by x*. Two-dimensional subspace orthogonal to 
the Killing vectors is described with p and z. The ranges of p and z are < p and 
— oo < z < oo respectively. The Latin indices i, j and k which run from 1 to — 2 
are used to assign the Killing vectors to the quantities {gij, ^ij, pi ...). The Latin 
indices p and q are kept for labeling the quantities with solitons. Both the indices 1 
and t mean time components. 

Especially in five dimensions, (j) and ip are used in place of x'^ and x^. And also 
in §3.1, the prolate-spheroidal coordinates x and y are introduced, which are defined 
by p = aV x"^ — — y^, z = axy with the ranges 1 < x and — 1 < y < 1- 

1.2. Rod structure 

The rod structure was first deliberately used to study the so called generalized 
Weyl solutions of D-dimensional vacuum Einstein equations, which admit D — 2 
orthogonal commuting Killing vector fields,'^ and next extended to the stationary 
and axisymmetric casesll^S In recent progress of solitonic solution-generations, the 
rod structure and practically the diagrammatic view of the rod structure become 
powerful tools. In fact, as examples shown in the following sections, using the analysis 
of rod structures we can guess and construct appropriate seed solutions to generate a 
desired solution. Also mathematical and physical properties of the resultant solution 
can be investigated from the viewpoint of rod structure with the aid of mathematical 
facts like the uniqueness theorems of solutions.'^ 

Originally 'rods' appeared in the process to construct the Weyl or generalized 
Weyl solutions as the rod sources of gravitational potential in hypothetical three- 
dimensional Newtonian gravity. For the static and axisymmetric case with the D — 2 
orthogonal Killing vector fields, the metric admits the diagonal form 

D-2 
1=2 

The functions Ui are axisymmetric solutions of the hypothetical three-dimensional 
Laplace equations 

{dl + Up + di^Ui = Q, (1-2) 

for i = 1, • • • ,D — 2, and also satisfy the constraint derived from the Einstein equa- 
tions 

D-2 

Y,U^ = \np. (1-3) 

i=l 

Once a set of {C/i, • • • , Ud-2} given, the quantity v can be determined by integrating 
the following equations 

^ 1=1 
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D-2 



(1-5) 



i=l 



From Eqs. ()l-2p and (jl-3p . the functions Ui can be considered as the Newtonian 
potentials produced by thin rods lying at some positions on z-axis, and all Ui must 
add up to the Newtonian potential corresponding to an infinite thin rod with linear 
mass density 1/2. We can say that when a set of rods satisfies the constraint p-3|) and 
each rod is assigned to one of Killing vector fields, a corresponding solution can be 
constructed uniquely. In this place, it should be, however, noticed that each rod does 
not necessarily have linear mass density 1/2, so that regularity of the corresponding 
generalized Weyl solution cannot be guaranteed. To construct a regular solution in 
which occurrence of conical singularities is allowed, each rod must also have linear 
mass density 1/2. When all of the rods have 1/2 density, equation ()l-3p means 
that there is only one rod corresponding to every points on the z-axis except some 
isolated points (i.e., junctions of adjacent rods). On each rod in this case, the metric 
coefficient e^' associated with the Killing vector field which is assigned to the rod 
behaves 0{p'^), so that the interval on the z-axis corresponding to the rod becomes 
a rotational axis or a horizon in the corresponding spacetime. 

We now introduce the rod structure of the solution as follows: first prepare the 
set of rods that was made by dividing z-axis into appropriate intervals; next assign 
a (D — 2)-dimensional vector (so called the direction vector) to each rod: 



The direction vector v in the generalized Weyl solutions just corresponds to one of 
Killing vectors On the interval corresponding to the rod, the norm of v is zero: 



Extension of rod structure to stationary and axisymmetric solutions with off-diagonal 
metric components can be done almost similarly: see Ref . I29p for complete discussion. 
In general case the direction vector v = v^-Mj can be defined through 



There being off-diagonal metric components, the direction becomes mixed with some 
Killing vectors. For example, the direction vector for the rod associated with a 
horizon is (1, ' ' ' , where Qi is a angular velocity of i-th rotation of the 

horizon, and for the rod associated with normal rotational axis the direction is 
(0, • • • ,0) where 1 appears in z-th slot (i = 2, • • • ,D — 2). Furthermore as 

an interesting but more complicated case, the rod structure of the spacetime with 
black lenses was also consider ed.'SS' 
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(1-7) 
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§2. Solution-generating technique I: Inverse scattering method 

2.1. Procedure of inverse scattering method 

In this section, we show the technique of constructing vacuum solutions, in- 
verse scattering method.'^''^ This technique can be applied only if we have D — 2 
commutable Killing vectors in D-dimensional spacetime. In the inverse scattering 
method, a known solution is transformed into a new solution. Here, we concentrate 
on how to construct vacuum solutions of spacetime and do not prove that vacuum 
Einstein equations are satisfied in the obtained new solutions. The detailed proof is 
shown in the original paper.® 

We begin with the D-dimensional stationary spacetime with D — 2 commutable 
Killing vectors (d/dx^), {d/dx"^), • • • , {d/dx^'"^), where (d/dx^) is a timelike Killing 
vector field and {d/dx"^), • • • , {d / dx^~'^) are spacelike Killing vector fields. Accord- 
ing to Earmark's paper,'2Sl' without loss of generality the metric can be written as 

ds^ = f{p, z){dp^ + dz^) + gij{p, z)dx^dx^, (2-1) 

where f{p,z) and gij(p,z) are a scalar function and an induced metric on (D — 2)- 
dimensional hypersurface, respectively, and the induced metric satisfies det gij = 
—p^. Einstein equations of this metric become 

dpU^' +d,V,^ = 0, (2-2) 
f^i' =p(5pto)/^ (2-3) 
^,'■ = p(5.to)/^ (2-4) 

d,lnf = j-U^'Vj\ (2-6) 

The first three equations ()2-2p - ()2-4p are a set of non-linear differential equations 
relating to only the induced metric gij. In the last two equations (j2-5|) and (|2-6|) . 
derivatives are operated on only the scalar function In /. In order to obtain solutions, 
we solve the first three non-linear equations for the induced metric gij, and then, 
integrate the last two equations. While the integration of the last two equations 
is not difficult, it is hard to solve the first three equations because of non-linearity. 
Applying the technique of solving the inverse problem in the scattering physics, we 
can solve the first three equations. 

In the inverse scattering method, we must prepare a seed metric which is an al- 
ready known solution of vacuum Einstein equations. Here, we suppose that we have 
a seed metric g^'^^ . The way to construct a seed metric is shown in the next subsec- 
tion. Defining a generating matrix p, z), we construct the linear differential 
equations in ^^j\x, p, z) as 
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First, we must solve these equations for since the new solutions are de- 

(1221) and (1231) arc 
than Einstein Eqs. (|2-2p - (|2-4p because they are linear equations in ^^^^ 



2Xp \ (0) _ ,^(0)^-Ay(f 
+ A2^^A j l^,, ^T^, (2-8) 

where U^^]^ and y ^'^^'^ are f/j and ^ made of gf^ , respectively, and A is a new 
complex parameter independent of p and z. 

XLU! 

scribed with and any parameters. Equations (|2-7p and (|2-8p are easier to solve 

1 Einstein Eqs. (1221)- ([231) becaus* 
Secondly, we introduce functions 

lip{p, z) = + {z- Up) - {z - ap), (2-9) 

fiqip, z)=-^p'^ + {z-aq)-{z- Uq), (2-10) 

where ap is a real constant. pp{p,z) and fiq{p,z) are called a soliton and an anti- 
soliton, respectively. In this section, we represent solitons pp and anti-solitons pp as 
p'p, collectively. We also introduce (D — 2)-dimensional vectors m^^^ associated with 
p'p. The (D — 2)-dimensional vectors m^P^ are called BZ vectors. 
Next, we construct n x n matrix as 

_ mf) i^'~'\p'p,p,z)-^fgl^ {¥^){p'q,p,z)-^fm^^ 

where n is the number of solitons and anti-solitons we introduce. 
Finally, a new metric is described as 

<7.=p-^"/(^-^)|^n^'f^-^^j^/',-, (2-12) 

9'^^ = (alf - E(^"')"V';V',-'ivf ^ivj"^) , (2-13) 

iv(^) = mf (2.14) 

The obtained new metric gij satisfies Eqs. ()2-2p - ()2-4p and det(7jj = —(?■ Actually, 
Eq. ([2431) also satisfies Eqs. (I22D-(l231). The factors P~'^'^l^^~'^^VCp=\P!]!^^~'^^ in 
Eq. ()2-12p are needed to satisfy det (7jj = —p^. 

A scalar function / can be also constructed with the scalar function of the seed 
metric /^*^^. Integrating Eqs. (|2-5p and (|2-6p with Eq. (j2-12p . we know the new 
scalar function / is 

( " 2(n+D-3) \ / " _4 A 

/ = C/(%i^ H/^'f n (/^'p-/^',)— detTp,. (2-15) 
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The inverse scattering method can not be apphed at the points where p = 
0. Equations (|2-2p - (|2-4p are singular on p = because of det^jj = —p^. The 
singularities are coordinate singularities or physical singularities. Thus, there is no 
guarantee that the new solution is regular on p = 0. After the transformation, the 
structures at points where p = must be checked. If we can find that there are 
no physical singularities at points where p = (or if we can remove all physical 
singularities by tuning the parameters), we can expected that the obtained metric is 
a regular solution of the vacuum Einstein equations. In the opposite point of view, 
we don't care about singularities at points where p = in a seed metric. In order 
to construct a regular metric, we pay attention to the structures only in the new 
metric. 

2.2. Diagonal seed metric 

In the inverse scattering method, we must prepare a seed metric. In this section, 
we show a manner of constructing seed metrics. Here, we restrict our attention to a 
diagonal seed metric because it can be easily constructed. 

We start from the following diagonal metric 

ffij^ = diag (-51,52,- •• ,9D^2)- (2-16) 

Then, each component is decoupled in Eqs. (|2-2|) - (|2-4p and we can solve them easily. 
The equations for the components are 

dpipdplngi) + d,ipdAngi) = (i = 1, ■ ■ ■ , D - 2). (2-17) 

These equations can be solved because they are linear equations in In^j. In fact, 
products of /o"", solitons and anti-solitons, that is, 

m 

gi = ±p^l[p'p, (2-18) 

p=i 

where n and m are arbitrary numbers, satisfy Eq. (|2-17p and hereinafter we use this 
solution. Because of det = — p^, diagonal metrics which are solutions of Einstein 
equations (except for on /? = 0) are represented as 

In sum, the way to construct diagonal metrics is to put minus sign (— ) at the first 
component, at a numerator of a component, and solitons (and/or anti-solitons) at 
numerators and denominators so that the number of the same soliton (anti-soliton) 
at numerators is equal to that at denominators. 

The scalar function of the seed metric f^^^ must be constructed. For ease in 
explanation of the method of constructing f^^\ we use the following notation 

/ ,,/("),/(")... ,/("),/(")... n2,/(") . 

(0) _ / l^ l,lt^ 1,2 M 2,lM 2,2 P M k,lh^ k,2 1 „„s 

% -diag - ' ,id) Ad) ' ,Ad),id) ^^■^^) 

\ A* l,l/^ 1,2 ■ ■ ■ 2,1^^ 2,2 ■ ■ ■ k,lf^ k,2' ' ' 
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where, of course, this metric is constructed in the way just described in Eq. ()2-19p 
and the A;-th component has the factor p^. Integrating Eqs. (j2-5p and (|2-6p with Eq. 
(|2-2U|) . we can obtain the concrete form of f^^\ 



P(0) 



,Ad) 
P k,p 



n 1 ( n '''' 

i=l [_ \p,q 



i,p i,q 




l(n) i(n) 

^ hP ^ hq 

p^q H ~T~ H- i^ph^ i^q 




( 



1^' 



in) 
i,p 



^ 2 , f /(") 



( 



n 



pjLqP +P'^,p^^^,q 

1 

2 



p + 



(2-21) 



where k is an integration constant. 



2.3. Inverse scattering method with diagonal seed metric 

Starting from a diagonal seed metric and performing a double transformation 
make the calculation of obtaining the new metric easier. Advantage of starting a 
diagonal seed metric is easiness to solve Eq. ()2-8p for 'I/ij. Advantage of performing 
the double transformation is that we don't need to care about the normalization 
factor in Eq. ()2-12p . We have an additional advantage that the seed metric is easily 
predicted. Due to these advantages, we often start from a diagonal seed metric and 
perform the double transformation. 

At first, we explain the double transformation. We denote the initial seed metric, 
the obtained metric after the first transformation and the final obtained metric as 
gfj\ (Jij^ and gij, respectively. In the double transformation we don't need to care 
about the normalization factor in Eq. ()2-12p . We will explain the reason later. 
We prepare some solitons /ip and anti-solitons fig and make anti-solitons p,p and 
solitons corresponding to the solitons fip and the anti-solitons Jig, respectively. 
In the first transformation, starting from a diagonal seed metric we use anti-solitons 
ftp and solitons with the trivial BZ vectors which are equal to (d/dt). The 
soliton transformation (without attentiveness to the determinant) with trivial BZ 
vectors is the same as an operation of a multiplication only in the tt component by 
U.pi-p'^/pDUqi-p'^/pl)- Thus, glf is written as 



~(o) 



diag Yl 



pI 



n 



-p \ (0) (0) (0) 
2 I 9tt '522 '533 
Pq 



(2-22) 



In the opposite point of view, the initial seed metric gf^^^ can be obtained by multi- 
plying gf^^ in only the tt component by 



Wi-pi/p") n(-^'/p^) = Wi-p'/pi) n(-pVA 



(2-23) 
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metric after 
first transformation 



(det 6 .7^ - p'-) 



new metric 

(det g = - p-) 



Fig. 1. Double transformation; Starting an initial seed metric glj \ we create the metric g^j' by the 
inverse scattering method with anti-solitons fip, solitons fiq and the trivial BZ vectors. Using a 
seed metric as the metric obtained in the operation of the first transformation g'^' , we perform 
the second transformation with solitons p.p, anti-solitons jlq and BZ vectors. Then, we obtain the 
new metric. If we use the trivial BZ vectors in the second transformation, the finally obtained 
metric is reverted to the initial seed metric. 



x(0) 



This means that the initial seed metric g\j can be obtained by the inverse scattering 

method with the sohtons fip and the anti-sohtons jlq starting from the metric g^j\ 
Therefore, we call the first transformation "removing the solitons fj,p and the anti- 
solitons fiq" . After that, we perform the second transformation. In the second 

transformation, starting from g^'j\ we use the solitons /ip, anti-solitons jlq and non- 
trivial BZ vectors. 

One of the advantages by starting from a diagonal metric is easiness to solve 
the linear differential equations (j2-7p and (j2-8p . In the inverse scattering method, 
we must solve the linear differential equations (|'2-7|) and (|2-8p for Because of 

the diagonalization of a seed metric, Eqs. ()2-7p and ()2-8p become diagonalized and 
can be solved easily. Complying the following guides, we can construct ^ij^ ■ We 

represent the seed metric S'^^H only with solitons and/or anti-solitons using the 
identity p'^ = —ppjlp and replace all solitons and anti-solitons fi'p by (/i'p — A). The 
obtained matrix satisfies Eq. (|2-8p . id est it is ^j^? ■ For example, when we deal with 



*' Since the first transformation is equivalent to the operation of the multiplication only in tt 
component, we don't need to create 5^'}°'. In the second transformation, the seed metric is the metric 

after the first transformation gl'^\ and thus, Eqs. p-7|l and H2-8|l is made of jji'j^ ■ We denote •f^'^'' 
as the solution of Eqs. (|2-7|l and (|2-8p in the second transformation when we concretely construct 
the solutions in §§3]and[S] 
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the seed metric (|2-20p . ^^j^ becomes 

■ (4"^-A)(41-A)--- (4";-A)(41-A) 



=diag| 



(4'} - A)(/.'g - A) • • • ' - A)(^'g - A) • • • ' 
-if^k,o - A)(Afc,o - A)(^tJ - A)(m11 - A) • • • 



(/^'S-A)(^'S-A) 



(2-24) 



Performing the double transformation removes the burden of normahzation in 
Eq. (I242H . The factors (^J]^^^ /i'2/(^-2)^ in Eq. (I242D are needed for 

tuning the determinant of a new metric gij to — /O^. This adjustment does not depend 
on BZ vectors, and thus, the BZ vectors do not affect the determinant of the metric. 
From the viewpoint of starting from the metric g^j^ , the initial seed metric gf^^ and 
the new metric gij are created by the transformation with the same solitons and 
anti-solitons while the BZ vectors are different. Thus, if we construct the initial 
seed metric g^'^^ whose determinant is — p^, that of the new metric gij automatically 
becomes —p"^- In the same reason, regarding the construction of the scalar function 
of the new metric / we do not need to mind the factors in Eq. (|2-15p and the new 
scalar function / is 

/ = C/W^^, (2-25) 
detr^°) 

where Fpg is created with gij\ the solitons /Xp, anti-solitons fig, and non-trivial BZ 

vectors, and Fpg^ is created with glj\ the solitons /Xp, anti-solitons p.g, and trivial 
BZ vectors. 

The initial seed metric g^^'^ belongs to the family of the new metric gij obtained 
by the double transformation. While the transformation with trivial BZ vectors by 
starting from the metric gf^^ gives the initial seed metric glj\ we can obtain the 
new metric gij by the transformation with non-trivial BZ vectors. Parameters of the 
transformation are introduced via the BZ vectors and, in this sense, the new metric 
gij becomes the initial seed metric g^'^^ by tuning these parameters. Therefore, it 

is expected that there are some similarities between the seed metric o-^^ and new 



metric gij. Actually, in the constructions of the solutions in §^ and [5l we start 
from the seed metric whose rod structure is similar to that of the new metric (but 
unfortunately not the same). 

2.4. Construction of metric corresponding to given rod structure 

Hereinafter, we concentrate on solutions in five-dimensional spacetime. More- 
over, since the inverse scattering method can be applied only to the case with three 
commutable Killing vectors, we suppose that we have them; (d/dt), {d/dcj)) and 
{d/dtp) where {d/d(p) and {d/dtp) are spacelike vectors and (d/dt) is a timelike 
vector. As explained in the last section, solutions of vacuum Einstein equations in 
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Fig. 2. Rod structure of example; linear mass density 1/2 is described as the thick black lines. The 
thick gray line means linear mass density —1/2. 



five-dimensional spacetime with three commutable Kihing vector fields can be repre- 
sented by rod structures. In this section, we explain how to construct the (diagonal) 
metric corresponding to a given rod structure. In a spacetime with a diagonal met- 
ric, the directions of the rods are (d/dt), {d/dcj)) or [d/dil)). Therefore, we suppose 
that all rods in the seed metric direct to {d/dt), {d/d(j)) and/or {d/dij)). In order 
to construct asymptotically flat spacetime with black rings (and a black hole), we 
prepare the seed metric only with linear mass density 1/2 and linear mass density 
—1/2. Thus, we consider the case where a seed metric has only such a rods. 

We explain the method to construct the seed metric from a given rod structure 
with an example described as Fig. [2l First, we write the minus sign (— ) in tt 
component. Next, we focus on the left end {z = —00) of the rod structure. We 
put in the numerator of the corresponding component to the existence of a rod 
at the left end. In the example, appears in cfxj) component. Moving to the 
right (increasing z) until the configurations of rods change, we construct the soliton 
corresponding to the changing point. If z = ai there, the constructed soliton is 

^^l = ^/p^ + {z-aif -{z-ai). (2-26) 

If a rod with linear mass density 1/2 ends on some component, we add the soliton 
in the denominator of the corresponding component. On the contrary, if a rod with 
linear mass density 1/2 starts, the soliton is added in the numerator. The soliton 
appears in the numerator of the corresponding component to the starting point of 
a rod with linear mass density —1/2, while it appears in the denominator of the 
corresponding component to the end point. Repeating this operation to reach the 
right end {z = 00), the obtained metric has the originally-provided rod structure. In 
the example, the obtained metric is written as 

Qij = diag , , . (2-27) 

V fJ-s--- f^2--- fJ-i--- J 

It is easy to make sure the obtained metric corresponds to the originally-provided 
rod structure. Rods with linear mass density 1/2 and rods with linear mass density 
— 1/2 mean the corresponding components become O(p^) and 0(p~^) in the p — )■ 
limit, respectively. In this limit, a soliton behaves 0{p^) for z < and 0{p'^) for 
z > ttp where ap is the constant parameter of the soliton. Adding the soliton pp in 
numerators is multiplying 0{p'^) only for z > ap, and thus, at the point z = ap the 
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corresponding rod with linear mass density 1/2 ends or the corresponding a rod with 
linear mass density —1/2 starts. Putting the soliton /Xp in denominators bring in the 
opposite effects. 

§3. Solution-generating technique II: Backlund transformation 

3.1. Backlund transformation 

In this section we briefly explain the solution-generating technique for the five- 
dimensional Einstein equations, which is called Backlund transformation. This is 
essentially the technique to generate new solutions of the Ernst equation from a 
known solution. By using this technique several singly spinning solutions for the 
five-dimensional spacetime were constructed .'^"'^^ '133 'EJ 

We start the analysis from the following form of the metric, 



(3-1) 

Using this metric form the Einstein equations are reduced to the following set of 
equations, 



(i) V^T = 0, 



f 5p7T = fp [{d,Tf-{d^Tf 



(ii) { 

(iii) V2£:5 

(iv) < 



V£s ■ V£s, 



dpls 
dzls 



P 



2(^5+^5) 



2(^:5+^5) 

(v) (9p^, d,^) = p-^e^^ i-dzco, dpuj) , 

(vi) 7 = 75. + 72., 

(vn) Ui = — , 



dpSsdpSs - dzSsdzSs ) 
dpSsdJs + dpSsdJs ) , 



where the function 'P{p, z) is defined through the equation (v) and the function Eg 
is defined by £3 := + i<P. It should be noted that and 'P corresponds to a 
gravitational potential and a twist potential. The equation (iii) is exactly the same 
as the Ernst equation in four dimensions. The most nontrivial task to obtain 
new metrics is to solve the equation (iii) because of its nonlinearity. We use the 
method similar to the Neugebauer's Backlund transformatiorl^ or the Hoenselaers- 
Kinnersley-Xanthopoulos transformation.'^^ ' 

Following the procedure given by Castejon-Amenedo and Manko,'23 for a static 
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seed solution a new Ernst potential can be written in the form 

5(0) x{l + ab) + iy{b — a) — (1 — ia)(l — ib) 
x{l + ab) + iy{b - a) + (1 - ia)(l - «6) ' 



where x and y are the prolate-spheroidal coordinates: p = aV x'^ — ly/l — y'^, z = 
axy with the ranges 1 < x and — 1 < y < 1, and the functions a and b satisfy the 
following simple first-order differential equations 



(lna),y 
(In 6),. 
(In 6),, 



—[{^y-i)sS^ + {i-y')sf], 

X — y ' 
^[-(x2-l)5W + (xy- 1)^(0)], 

(xy + 1)5(0) + (1-2/2)^(0) 
(x2-l)S(0) + (x2/ + 1)5(0) 



x-y 
1 

x + y I 
1 



x + y 



(3-2) 



The corresponding expressions for the metric functions can be obtained by using the 
formulas shown by Ref. [37j) . For the seed, 



+ e 



(3-3) 



a new solution is given by 

S'(o)_T(0) ^ 



B 



dt-{2cje-''''^ + C, 



, -sW-TC) 2/ 2 2n 7,2 

+— e cr (x - 1)(1 - y )d4) 



'sm-rm ^^2^'^ -V^ ( dx 



+ 



dy' 



x-^ — 1 \x^ — 1 1 — 
where Ci and C2 are constants and j4, B and C are defined by 
A := (x^ - 1)(1 + abf - (1 - y2)(& _ af , 



(3-4) 



C7 



[(x + 1) + (x - l)ahf + [(1 + y)a + (1 - y)6]2, 
(x2-l)(l + G6)[(l-y)6-(l + y)a] 
+(1 - y2)(6 _ a)[(x + 1) - (x - l)ab]. 



(3-5) 



The function 7' in Eq. (|3-4|) is a 7 function corresponding to the static metric, 



, 2 -T(0) 
as = e 



+e 



2(y-C/r')-5(0) (J^2^^^2y 



+ e 



(3-6) 
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where f/n^^' = — In 



1 



x + 1 



Next we consider the solutions of the differential equations (|3-'2p . At first, we 
examine the case of a typical seed function 

S(^^ = lln[Rd + {z-d)], (3-7) 



where Rd = y^p^ + (z — d)^. The general seed function is composed of seed functions 
of this form. Note that the above function S^'^^ is a Newtonian potential whose source 
is a semi-infinite thin rodPSJ 

In this case we can confirm that the following a and b satisfy the differential 
equations (j3-2p . 



where 



a = I, 



Od,c 



2 



(3-8) 



(3-9) 



Here the functions Ud and Uc are defined as Ud '■= \ [^d + (-^ — c^) ] and Uc '■= 
^ In [ — (z — c) ] . Because of the linearity of the differential equations (|3-2p for 
we can easily obtain a and 6 which correspond to a general seed function if it 
is a linear combination of ()3-7p . 

The function 7' is defined from the static metric (|3-6p . so that 7' obeys the 
following equations, 

dpi' = \P [{dpS'f - [d.S'f] + \p [{d.T'f - [d.T'f] , (3-10) 

d.i = \p {d,S'd,S'\ + \p {d,T'd,T'\ , (3-11) 

where the first terms are contributions from Eq. (iv) and the second terms come 
from Eq. (ii). Here the functions S' and T' can be read out from Eq. (j3-6p as 



5' = 2C/f^)+5W, 



(3-12) 
(3-13) 



To integrate these equations we can use the following fact that the partial differential 
equations 



dpl'cd = P dpUcdpUd - d^UcdJJd 



dzicd = P 
have the following solution, 

led = 



dpUcdzUd + dpUddzUc 



^Uc + ^Ud-\lnYcd, 



(3-14) 
(3-15) 

(3-16) 



where Y^d ■= RcRd + {z — c){z — d) + p'^ . The general solution of 7' is given by the 
linear combination of the functions 7^^. 
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3.2. Relation between two solution-generating methods 

To construct five-dimensional solutions, we have been mainly using two differ- 
ent solitonic solution-generating techniques, the inverse scattering method and the 
Backlund transformation. Both methods successfully constructed several important 
solutions. Also, these two methods will be used to find new higher-dimensional black 
hole solutions. Note that while the solutions with single angular momentum com- 
ponent can be constructed by the both methods, the solutions with double angular 
momentum components can be constructed only by the inverse scattering method. 

It is important to consider the relation between solutions generated by these two 
methods. The relation between the four-dimensional solutions generated by these two 
methods was investigated in Refs. [T]) J38 |) -H0 |) . Recently the relation between the five- 
dimensional solutions was also investigated by Tomizawa et al^^ The singly spinning 
five-dimensional solutions found or reconstructed by the Backlund transformation 
can also be constructed by the inverse scattering methods. The solutions obtained 
by Backlund transformation are 2-soliton solutions whose seed metrics are diagonal. 
It was shown that the 2-soliton solutions generated by the inverse scattering method 
coincide with ones generated by the Backlund transformation for the general diagonal 
seed. In fact, the five-dimensional black ring with rotation was constructed from 
the Minkowski spacetime by using both methods .KJ Also it was shown that the 
black ring with rotation can be constructed from the Euclidean C-metric by using 
both methods.13'1^ 

It has been often used the inverse scattering method that was modified by Pomer- 
ansky to construct the five-dimensional solutions.'^ In single-rotational case we can 
easily transform the procedure of the original inverse scattering method into the 
modified method. This means that we can investigate the relation between solu- 
tions obtained by the Backlund transformation and the modified inverse scattering 
method through the mediation of the original inverse scattering method. 

§4. Single horizon solution 

Previously five-dimensional black hole solutions which have only one event hori- 
zon were found without the use of solitonic solution-generating techniques. These 
solutions were reconstructed and generalized by using the solitonic techniques. In 
this section we briefly explain the ways to construct the solutions of five-dimensional 
Myers-Perry black hole and black rings. 

4.1. Myers-Perry black hole I 

The Myers-Perry black hole solution with single angular momentum component 
was regenerated from the Minkowski seed by using the solitonic techniques. This 
solution is obtained as a limit solution of S"^ rotating black ring without the inside 
space of the ring.im'''I3' Also it can be directly derived from the same seed of S'^ 
rotating black ring by adding solitons at appropriate positions. 

The doubly spinning solution of Myers- Perry black hole was regenerated by using 
the inverse scattering method."^ To obtain the seed metric of the solution we start 



t 



— oo (l\ 0,2 00 

Fig. 3. Rod structure for the seed metric p'"-' of Myers-Perry black hole I. Solid lines corresponding 
to sources of linear mass density 1/2. 

from the Schwarzschild-Tangherlini black hole 

^7j?=diag(-^/,^,) (4-1) 

and 

The rod structure of this metric is shown in Fig. [3j We divide the tt component of 

2 2 

gio) \)y _E:2 i^y which we effectively remove an anti-soliton at 2: = ai and a 

p Pi 

soliton at z = 02 from g^'^\ We denote the obtained seed metric as g' , 

^2 



9f=di-4-^A,,.]. (4-3) 
For the convenience, we use the seed metric which is rescaled by a factor of — — , 

=diag(l,/22,-m)- (4-4) 

The corresponding generating matrix ^^^^ which satisfies Eqs. (j2-7p and (j2-8p is 
given by 

= diag (1, (/I2 - A), -(/ui - A)) . (4-5) 
We perform a 2-soliton transformation with g^^^ where an anti-soliton is placed at 
z = ai with BZ vector ttiq^^ = (1,6,0) and a soliton at z = 02 with = (1,0, c). 
We denote the resulting metric as g and rescale it to find the metric g = —j^9- 
There is a freedom to make linear transformation in the space of coordinates t, 4> 
and ip. For example, we perform a following linear coordinate transformation 

t = + 26(02 - ai)-/-"''^ + 2c(a2 - ai)V'°''™ (4-6) 
(j) = 2(ai - a2W^'" + 6cV'°''™ (4-7) 
iP = 2(ai - a2W^'" + 6c0°''" (4-8) 

to remove global rotations of the solution. The conformal factor of two-dimensional 
part / is given by 



where F is found in the process of constructing g and F^'^^ = F\ 



b=c=0- 
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Fig. 4. Rod structure for the seed metric g'''' of Myers- Perry black hole II. Black lines corresponding 
to sources of linear mass density 1/2 and gray lines to linear mass density —1/2. 



4.2. Myers-Perry black hole II 

There is another way to construct the doubly spinning solution of Myers-Perry 
black hole by using the inverse scattering method. We start from the following seed 
metric 



and 



g]-' = diag , , 4-10 

P4. /U2/i3 Pi J 



r{0) _ P2P3 W12W13WUW24W34 ,^ 

^ ~ Pi WIWUW22W33W44 ^ ■ ' 

where Wpq = + PpPq- The rod structure of this metric is shown in Fig. [H At first 
we remove a soliton at 2; = ai and an anti-soliton at 2; = 02 with trivial BZ vectors 
(1,0,0) as 

-■1?=diagf-^,^,^V (4-12) 



Pi P2P3 Pi J 



For the convenience, the seed metric is rescaled by a factor of 



tJ.4 ' 



-(°)-diag(l,^,-^). (4-13) 
P3 PA J 



After construction of the generating matrix as 

^(0) ^ di^g ( , {-P2 - X){t.l - X) ^ _ {p2-X)iP3-X) \ 

\ P3- X Pa- X J 

we perform a 2-soliton transformation with g^^^ where an anti-soliton is placed at 
z = ai with BZ vector mg^^ = (1, 0, c) and a soliton at z = 04 with TOq^^ = (1, b, 0). 



Finally we rescale the metric g to find the metric g = —^9- The conformal factor 
/ is given by 



where F^^^ = r\h=c=o- The BZ parameters should be determined to remove the 
singularities on the 



_^ / 2(a4 - ai)(a4 - 02) ^ _ _^ / 2(a2 - ai)(a3 - Qi) (^^.IQ) 
V («4 - 03) ' y (04 - ai) 
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Fig. 5. Rod structure for the seed metric g'"' of S'^ rotating black ring. Solid lines corresponding 
to sources of linear mass density 1/2. 



The obtained metric {g, f) already has the same asymptotic form as the Minkowski 
spacetime. Therefore we do not need to perform a linear transformation. 

4.3. Black ring 

The five-dimensional black ring has an event horizon whose topology is x S"^. 
The black ring has to rotate along the direction to avoid a conical singularity 
in the inside space of the ring. The original solution of black ring which has only 
one angular momentum component along the direction was found without use 
of solitonic techniques.® The black ring solutions which have angular momentum 
component along the S"^ direction were constructed by using the solitonic solution- 
generating techniques. 

4.3.1. Black ring with S"^ rotation 

The solution of black ring with S'^ rotation was constructed by using the Backlund 
transformatioiPD''^ and reconstructed by using the inverse scattering method.'^ 
Also the same solution was found by the educational guess work.SSJ The seed metric 
of S'^ rotating black ring is the five-dimensional Minkowski spacetime when we use 
the Backlund transformation. While we can construct the same solution from the 
same Minkowski seed by using the inverse scattering method, it is convenient to start 
from the metric of static black ring, 

,(0) _ „ f W 



5.7=diag(-^,^,/x3) (4-17) 

and 



= ^3 "^^f^' . (4-18) 

The rod structure of this seed metric is shown in Fig. [5j At first we remove an anti- 
soliton at z = ai and a soliton at ^ = a2 with trivial BZ vectors and next rescale it 
by a factor —j^- The resulting metric is 

5;?=diag(l,/l3,-^). (4-19) 
After constructing a generating matrix from this metric, 

= diag (l, - A, , (4-20) 
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we perform a 2-soliton transformation to obtain the metric g by adding an anti- 
soliton sX z = ai with a BZ vector (1, 61, 0) and a sohton sX z = a2 with a BZ vector 



(1, 62) 0) to the seed Q^ij^ ■ Finahy we rescale the metric g to find the metric g = —^9- 
The conformal factor / is given by 



4(ai-a3)^ (0) detr 
(2(ai -as) -6162)2^ detrW ^ ' 

where F^^^ = -r|f,^=fe2=o- In order for the metric to be asymptotically Minkowski, we 
need perform an appropriate linear transformation of the coordinates t and 

t = t"""^ - ^("l- "2)^l ^ncw^ ^ ^ ^new_ .4.22) 

2(ai - as) - 6162 

In general the solutions has some undesiered features. When the parameters satisfy 
the following equation 

2(a2 - as)6i + 2(ai - as)62 - bib2{bi + 62) = 0, (4-23) 

the finite spacelike rod between z = a2 and z = as corresponds with the (j) axis 
where g^ptj) = 0. Even in this case, we can not remove conical singularities in the 
inside space of the ring. 

4.3.2. Black ring with rotation 

The black ring with an angular momentum component along direction was re- 
constructed by the Backlund transformatiorf^ and the inverse scattering methocP^ 
starting from a Euclidean C-metric as a seed. In these analyses, a 2-soliton transfor- 
mation was used to derive the solution. We can also construct the rotating black 
ring by a 1-soliton transformation from another seed^^^^ To do so, we start from 
the following diagonal metric, 

9ij = diag , , . 4-24 

The conformal factor of the seed is given by 

.(0) ^ /^2^4^f2^13W-15W-i6l^23^25^|4^36W|5^56 .^.25) 
m WuW24Wi^Wi^W^eULlWii 

The rod structure of this metric is shown in Fig. [6l We need a negative density rod 
between z = ai and a2 to obtain the rotating black ring. We remove an anti- 
soliton at z = ai with trivial BZ vector and multiply the factor ^ to the metric. 

The resulting metric is denoted as g^^^ 



gl-' = diag — , , -— . 4-26) 

Then we construct a generating matrix from g^'^^ as 

^(■^^diagf^. fe-;''"'-;;' -Igi^l, (4.27) 

\fi3-X (^1 - A)(//4 - A) /X4-A 
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Fig. 6. Rod structure for the seed metric i;'"' of rotating black ring. Black lines corresponding 
to sources of linear mass density 1/2 and gray lines to linear mass density —1/2. 



Next we readd an anti-soliton with a nontrivial BZ vector (1,0, c) and obtain a metric 

2 

g. Finally we rescale the metric g to obtain the metric g = ^9- The conformal factor 
/ is given by 

where T^o) = r|c=o. To remove the singular behavior on the axis we set the BZ 
parameter as 

' ^^K^^^. (4.29) 




04 — ai 

In this case we do not need to perform a linear transformation of the coordinates. In 
general the solution has a conical singularity on the plane of the ring. The balance 
condition to avoid a conical singularity is given by, for example, 

al - 2a2ai + 0304 

ai = — . (4-30) 

03 - 04 

4.3.3. Doubly spinning black ring 

The doubly spinning black ring solution was found by using the inverse scatter- 
ing method.!^ The seed metric of this solution was obtained by removing a pair 
of solitons with trivial BZ vectors from the Emparan-Reall black ring which is ro- 
tating along direction. Therefore the seed metric is not diagonal. To rotate the 
Emparan-Real black ring along the S"^ direction we readd the same pair of solitons 
with nontrivial BZ vectors to the nondiagonal seed. The most disturbing difficulty 
is to find a generating matrix field iP''^''^ which corresponds to the nondiagonal met- 
ric of Emparan-Reall black ring. Fortunately, it is known that if we construct the 
rotating black ring via a one-soliton transformation as in 14.3.21 this matrix field 
automatically is given. It is worth to note that the regular black ring with two an- 
gular momenta comes from the seed solution, which is obtained from the regular S"^ 
rotating black ring without conical singularity. 

§5. Solution of disconnected horizons 



In five dimensions, in addition to the solutions which have only one event hori- 
zon component, there exist solutions which have more than one disconnected event 
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Fig. 7. Rod structure for the seed metric g'^^^ of black Saturn. Black lines corresponding to sources 
of linear mass density 1/2 and gray lines to linear mass density —1/2. 



horizons. These solutions were obtained by using solitonic techniques, the Backlund 
transformation and the inverse scattering method. In the fohowing subsections we 
win briefly explain the ways to construct these solutions and the physical features 
of them. 

5.1. Black Saturn 

One of the example of black hole solutions with disconnected horizons is the 
black Saturn solutionj^^ in which a spherical black hole is surrounded by a black 
ring. This solution was constructed by using the inverse scattering method starting 
from the following seed metric 



no 

and 



.S^=diag(-^,^,^) (5.1) 



The rod structure of this metric is shown in Fig. [71 At first we remove anti-solitons 
at z = ai and 04 and a soliton at z = 05 with trivial BZ vectors and rescale it by a 
factor of to flnd 

^r=diag(J-,^,-^). (5-3) 
The generating matrix is given by 

^(0) = diag ( ' . (5.4) 

V/U3 - A (/i2 - A)(/X5 - A) /U2 - A y 

Next we perform a 3-soliton transformation with 5(0) as seed: Add an anti-soliton 
at 2; = oi with BZ vector (1,0, ci) and an anti-soliton at z = 04 with (1,6,0) and a 
soliton at z = 05 with (1, 0, C2). We denote the resulting metric g. Finally we rescale 

2 

q to find a metric q = ^^-^q. The metric factor f is obtained as 



where -T^"^ = i^l 



6=ci=C2=0- 
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Fig. 8. Rod structure for the seed metric black di-ring which is used in the analysis by Evslin and 
Krishnan. Black lines corresponding to sources of linear mass density 1/2 and gray lines to 
linear mass density —1/2. 



The physical feature of black Saturn was investigated only for the case of single 
spin solution in which 6 = 0. One of the most important properties of multi- horizon 
solutions is continuous non-uniqueness of the solution. To show this property, the 
phase diagram of the black Saturn was investigated in Refs. ri8|l . Hl|) . The plot of 
random sets of points in the phase diagram showed that the black Saturn covers the 
wide region of the phase diagram. This results can be understood by the analysis 
based on the thin and long ring approximation in which the black Saturn can be 
modeled as a simple superposition of an Myers-Perry black hole and a very thin 
black ring.SS It was argued that the configurations that approach maximal entropy 
for fixed mass and angular momentum are black Saturns with a nearly static black 
hole and a very thin black ring. 

Another important property of multi-horizon solutions is the existence of ther- 
modynamical equilibrium configuration in which the horizons have the same temper- 
ature and angular velocity. Requiring equal temperatures and angular velocities 
for the central black hole and the black ring imposes two conditions on the parame- 
ters. This removes entirely the continuous non- uniqueness, leaving at most discrete 
degeneracies. Black Saturns in thermodynamical equilibrium thus form a curve in 
the phase diagram. The black Saturn in thermodynamical equilibrium is not ther- 
modynamically stable because the entropy of it is not the maximum. It has been 
shown that metastablity occurs when the dimensionless total angular momentum 
lies in a narrow window 0.85483 < < 0.85494 of the thin ring branch 

5.2. Black di-ring 

5.2.1. Construction of black di-ring 

We can also construct a multi-ring configuration by using solitonic solution- 
generating techniques. The first example of multi-ring solution is a black di-ring 
which has two concentric black rings rotating on the same plane. The black di- 
ring solutions were first constructed by using the Backlund transformatiorf^ and 
reconstructed by the inverse scattering method; The representations of these 
two solutions are very different. One of the reason for this difference is that these 
two representations are obtained by solitonic techniques starting from different seed 
metrics. 

In this subsection we derive a representation of black di-ring by using the inverse 
scattering method starting from a seed metric which is different from the one used in 
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Fig. 9. Rod structure for the seed metric g'*^' of black di-ring. Black lines corresponding to sources 
of linear mass density 1/2 and gray lines to linear mass density —1/2. 



the analysis of Evslin and Krishnan. This representation of black di-ring directly cor- 
responds to the one obtained by Backlund transformation because the seed metrics 
are essentially same. The representation obtained by Evslin and Krishnan is pre- 
sented in Chap. 2, which is derived from the seed metric whose rod structure is given 
in Fig. [8l Considering the relations between these solutions, it was shown that the 
different solution sets of di-rings which are obtained by the Backlund transformation 
and the inverse scattering method are completely equivalent 

The seed solution is represented by the rod structure given in Fig. [9l The 
corresponding seed metric is given by 

= diag (-^, (5-6) 



and 



(0) ^ fcV2^4^7 Wl2Wl^Wl^WnW2^'W2&W^AW^^W:i^^W^^W^f,''WQ^ 

The solution of black di-ring is constructed as follows: 

1. Perform two 1-soliton transformations on the seed solution: Remove an anti- 
soliton at 2; = fli with trivial BZ vector (1,0,0) and a soliton at 2; = 04 with 
trivial BZ vector (1,0,0). The result is the metric matrix 

= diag f-^, . (5-8) 

2. Rescale g'^^'^ by a factor of — — to find 

/i4 

^W=diagf^,^^i^^,-^V (5-9) 
V/U6 fJ'ifJ-bfJ-r ^J'3 ) 



2 

where = — — • This will be the seed for the next solution transformation. 
3. The generating matrix can be found from It is 



•^(°)(A,p,z) = diag 



(/X5 - A) (/ii - A)(//2 - A)(;U3 - A)(/i6 - A) 

(5-10) 



(/X6 - A) ' (//4 - A)(/Li5 - A)(//7 - A) 

(//2 - A) (^7 - A )^ 

(/^3 - A) 











(1,02,0) 
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Fig. 10. Rod structure of the regular black di-ring. Direction vectors are presented over each rod. 
Singularities at 2: = ai and 2 = 04 are removed by setting the BZ parameters ci and C2 correctly. 



4. Perform a 2-soliton transformation with g^'^^ as a seed: Add an anti-soliton at 
z = ai with BZ vector ttiq^^ = (1, 0, ci) and a sohton at 2: = 04 with BZ vector 

ttiq = (1,0,02). Denote the resulting metric g. 

5. Rescale g to find 

5 = -^5. (5-11) 

6. The conformal factor / is constructed using Eq. (j2-25p . 

The result ((7, /) is the final solution after performing an appropriate linear coordi- 
nate transformation. The exact expressions of the final metric are written down in 
the next subsection 15.2.21 The rod structure of regular black di-ring with single spin 
is given in Fig. [TOl 

5.2.2. Solution of black di-ring 

We write down the metric functions of black di-ring. The non-zero components 
of the metric can be written as 

^ Xi + ci^Xa + C2^X3 + C1C2X4 + ci^ca^Xs 

{Di+Ci^D2 + C2^Ds + CiC2Di + Ci^C2^D5y ^ ' 



Zi + Ci^Z2 + C2^Z3 + C1C2Z4 + Ci^C2^Z5 
m^l-S {Di + Ci'^D2 + C2'^D3 + C1C2D4 + Ci'^C2^D^) ' 



(5-15) 



/U2/i5At7 



and 



. ^ C2^ i^l + Ci^£>2 + C2^£'3 + C1C2-P4 + Cl^C2^Z^5 .(Q) 

^ (Cl-C2)2 /,|/Z4;uiM|4M|6M|7PFi42PFi22iyi62iyi^ ' ^' ^ 



The definitions of function L), 's are 



= ^4^lA^ilM2A^M^(,''M„''Wl2^Wl^Wl(s^Wn'', 



(5-18) 
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L>2 = pVlA^2M3M4A^5M6M7Ml42M24^M462M47^VFl3^VFl5^ (5-19) 

D3 = pV2/^3/X4/X5M6Ai7Ml4^M34^M45^^yl2^H^16^H^17^ (5-20) 
D4 = 2/Xi/X2/i3/^4/^5/^6/^7M24M34M45M46M47l^llT^12T^13W44H^15W^16W^17, (5-21) 

and 

D5 = fill4l^4f4fJ'7M3A^M45^Wn^Wu''Wi5^ (5-22) 

where 

Mpq = fJ-p- fJ-q. (5-23) 

The functions Xj's are defined by the function DiS as 

X, = -t^D^, (5-24) 

X2 = ^D2, (5-25) 

/il/i4 

X3 = ^Ds, (5-26) 

X4 = -AI5-D4, (5-27) 

and 

X5 = -^L>5- (5-28) 
Ml 

The functions i^'s are defined as 

= -fi2f^3fi4f4f^7MuM24'M^e^Mir^WnWi2Wi3WuWi5Wi(iWi'r, (5-29) 

>2 = /X2/X3/i5/^7Mi4M24M34M45M46M47W^12^W^14W44VFi6^VFi7^, (5-30) 

^3 = -/X?/xi^5/^6/X?Ml4M34M24M45M46M47VFl3^W^14W44VFl5^ (5-31) 

and 

Y4 = /xi/X4/^5M6/i7^14M34^M45 ViiM^12Wl3Wi4W^15Wi6VFl7. (5-32) 

The functions ZiS are defined by the function Dj's as 

Zi= ji2iiAH7Di, (5-33) 

Z2 = -M^^,, (5.34) 

^3 = -^i^3, (5-35) 

^4 = /^l/^2/^7-D4, (5-36) 

and 

= ^^!^D,. (5-37) 

/X4 
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To remove singularities at z = oi and z = 04, we set the BZ parameters 



a3i«5i 



/ 2042064074 „„x 

C2 = ±W , (5-39) 

V 043054 

where apq = Op — Og. 

Note that, in order for the metric to asymptotically approach the Minkowski 
spacetime without global rotation, we have to replace and t(p components of the 
metric by using the freedom of linear coordinate transformation as, 

g^xj,-^ 9<l>4> + '^Cgt^ + C^gtt, 9t<t> ^ 9t4> + C gu (5-40) 

where 

C — 204104304502 (5-41) 

2042 064074 — 0430540102 

Using the exact expressions of black di-ring, we can compute physical variables of 
black di-ring. For example, the angular velocities of outer and inner horizons which 
appear in the direction vectors of rod structure of Fig. [TD] are obtained as 

= («42C1+021C2)(C1-C2) 

04i(2(042Ci + O21C2) - CiC2(Ci - C2)) 

and 

fi. = -r, , , (5.43) 

2 ((a64a74Cl - 06l07iC2)(Ci - C2) + a|^CiC2j 

5.2.3. Physical properties of black di-ring 

In this subsection we briefly introduce the results of investigations for the phys- 
ical properties of black di-ring. As similar as the black Saturn, one of the most 
important feature of black di-ring is continuous non-uniqueness. This fact can be 
confirmed by a plot of random systematical sampling in the phase diagram. 

Previously, it was argued that it is unlikely that multi-ring solutions exist as 
states in thermodynamical equilibrium because a black ring is uniquely determined 
by fixing the temperature T and angular velocity 17. Contrary to this expecta- 
tion, it was shown that the thermodynamical equilibrium black di-ring is possible.l^S 
Requiring equal temperatures and angular velocities for the both horizons imposes 
two conditions on the parameters. As a result, thermodynamical equilibrium black 
di-ring form a curve in the phase diagram. 

The black di-ring in thermodynamic equilibrium is not thermodynamically stable 
because the entropy of it is not the maximum. In addition, a possibility that the 
thermodynamic equilibrium is local maximum of the entropy is denied by the numeric 
search.!® While the window of the metastability is closed, it was found that there 
IS a narrow window 0.92075 < f < 0.92084 of the thin ring branch where the both 
eigenvalues of Hessian matrix of the entropy function are positive. In this region the 
entropy of the thermodynamic equilibria is locally minimum. It is likely that the 
interaction between the black rings makes the system thermodynamically unstable. 
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Fig. 11. Rod structure of seed metric of orthogonal black di-ring: rods with linear mass density 
1/2 are described as the thick black lines. The thick gray hue means linear mass density —1/2. 



5.3. Orthogonal black di-ring solution 

Here, we construct the solution with two black rings whose planes of rotations 
are orthogonal to each other. In the five-dimensional spacetime, we can take two 
spacelike orthogonal planes. In the orthogonal black di-ring solution, on each plane 
a rotating black ring exists. 

The orthogonal black di-ring solution was constructed with the seed solutiorP^ 

yj? = diag (-^, i^^] (5.44) 

and 

, /i2/i4/U6 Wi2Wi3Wi4WiQWi7W23W25W27W34W35W3eW4,5W47W5(iW57We7 



i^m> w,,^w,,^w,,^w,,^ul=iWpp 

(5-45) 



whose rod structure is described as Fig. [TTl Removing a soliton ^7 and an anti- 
soliton fli with trivial BZ vectors and rescaling it by a factor ^1/^7, we have 

ffW = diag f-^, . (5-46) 

and the generating matrix is written as 



<^(°) = diag 



(/»5 - A) (/Ui - A)(/i3 - A)(/Z4 - A) 
(/"3-A)' (/i2 - A)(/i6 - A) 

(/X2 - \){lJ.Q - X){jl7 - A) 
(^5 - A)(/i4 - A) 



(5-47) 



Then, we perform a 2-soliton transformation with 5(0) as seed metric. Here, we 
use a soliton ^7 with a BZ vector (1,6,0) and an anti-soliton Jii with a BZ vector 
(1,0, c). Rescaling the obtained metric by a factor ^17/^11, we have the orthogonal 
black di-ring solution. The scalar function / is obtained as 



where -T^"^ = i^l 



b=0,c=0- 
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The metric of the orthogonal black di-ring solution is written as 



/ = hH/F, (5-49) 
H = F + b^F^,) + c2F(,) + b'c'F^^,) , (5-50) 



gtt = H-\A + + cM(,) + 6^cM(6e)), (5-55) 

/9Vi/^7M572W-152^nVF77 ^ ^ 

^''^ pV2/"3Ai6M372vy,52pF472VFilW^77' ^ ' 



4^ ^ _ A^1A^7-'"14 -'"27 -'"67 ''''15 /tr crQ\ 



--B + b'H-\B^,)+c'B^,,-)), (5-60) 
S = ^^, (5-61) 

Ai2A*4^6 

pVl^3;^5M272Mg72M/,32 

""(^^ /.i/.iA^.7M3,2H/,^2^,,2^n' 

g^^ = C + c^H-\C(^,) + 62C(fee)), (5-64) 
C=^^^, (5-65) 

_ ^i^i^i^7Mi42M372W-i52 
^ ^1^2/X3/X4/X6/^7Mi42M272Mg72Vt/^^2 
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gt^ = H-'h{D + c'D(^,)), (5-68) 
gt^ = H-^ciE + b^E(,)), (5-71) 

This solution has generally singularities at (p, z) = (0, oi) and (0, ay) and conical 
singularities on the rods z G [—00, 02], G [osi 04], z G [04, 05] and z € [oe, 00]. These 
singularities can be removed if we set the BZ parameters 



'2071073^074 

a72«75«76 



(5-75) 



/ 2021031061071 tr^n\ 

and the integration constant in / 

= 1, (5-77) 

and if we impose the balance conditions 

_ 014016017025027034035036 , s 

i — 2 — 2 — 2 — 2 ' 

"15 '^24 ^26 ^37 

074072071053061054053052 r7n\ 

1 = 2 2 2 2 • (^'^^j 

O73 Og4 Ogg Og-i^ 

The rod structure of the obtained regular metric is illustrated in Fig. [121 The 
directions of rods on horizons are 

(l,^2ll\^2^^'>^ on zG [02,03], (5-80) 
fl,r?i'\r?f) on zG [04,05], (5-81) 



where 



.(1) _ 


b 075 076 

2 Oyg^ O7I 




,(2) _ 


6075 

2 071 074 ' 





= ^ _^^^5^ (5.82) 

^ 031 O61 071 



ZO61 071 





(i.nw.nw) , ( 
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Fig. 12. Rod structure of regular orthogonal black di-ring: rods with linear mass density 1/2 are 
described as the thick black hnes. 



This regular solution has four parameters which are corresponding to the radii and 
the angular momenta of the rotations of both black rings. Each black ring rotates 
to the S"^ direction because of the drag by the rotation of the other black ring. 

5.4. Black multi-ring 

It is expected that we can construct the solution where there are superposed 
multiple rotating black rings. Here, we show the expected method to construct 
the solution. 

We try to construct the solution where m + n black rings exist. We consider 
the situation where planes of m black rings are the same and orthogonal to the 
planes of the other n black rings. From the construction of rotating solutions 
in the previous subsections, it is suspected that we should introduce the non-trivial 
structure outside of each black ring in the seed metric. Then, the seed metric is 
described as Fig. [13] and written as 



p 



m n 9 

1 /^P3 l^ql 1^0 ^J-^ /^p2 ^J^ Atg3 ■ f^pl Atg2 



=diag -n— n— '-n— n— '/^oii— n— 



where fiqi is the soliton constructed with the constant Ogj. 

Starting from the seed metric Eq. ()5-84p . we remove m anti-solitons fipi {p = 
1 ■ ■ ■ m) and n solitons /igi {q = 1 • • • n). After that we perform the transformation 
with m anti-solitons p,pi and n solitons /igi where BZ vectors associated with the 
anti-solitons fipi and with the solitons Hqiave (1,0, Cp) and (1,6^,0), respectively. 
Generally, the obtained solution has singularities at z = Opi and z = a^i and conical 
singularities on the axes. It is expected that these singularities can be removed by 
tuning the parameters bq, Cp, Opi and a^i and that the regular solution has 2m + 2n 
parameters. They are corresponding to the radii and the angular momenta of the 
5^ direction of m -|- n black rings. 
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Fig. 13. Rod structure of seed metric of black multi-ring: rods with linear mass density 1/2 are 
described as the thick black lines. The thick gray line means linear mass density —1/2. 
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